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How much can we know about our Universe? All of our observations are restricted to a finite 
volume, and therefore our estimates of presumably global cosmological parameters are necessarily 
based on incomplete information. Even assuming that the Standard Model of cosmology is correct, 
this means that some cosmological questions may be unanswerable. For example, is the curvature 
parameter Q.k positive, negative, or identically zero? If its magnitude is sufficiently small, then 
due to cosmic variance no causal observation can ever answer that question. In this article, we first 
describe the gauge problems associated with defining the cosmic variance of cosmological parameters, 
then describe a solution involving the use of parameters defined on the surface of last scattering, and 
finally calculate the statistical variance of ideal measurements of the matter, radiation, and curvature 
density parameters. We find that Q.k cannot be measured to better than about 1.5 x 10^"" (Icr), and 
that this limit has already begun to decrease due to the flattening effect of dark energy. Proposed 
21 cm hydrogen experiments, for example, make this limit more than just a theoretical curiosity. 

PACS numbers: 98.80.Jk 



I. INTRODUCTION 

The Standard Model of Cosmology is a remarkable 
scientific achievement. The combination of general rel- 
ativity with statistical homogeneity and isotropy yields a 
model which currently requires only about a half dozen 
parameters [1] in order to provide a successful description 
of the real Universe. 

Assuming this Standard Model is correct, the physical 
Universe is larger (possibly infinitely so) than the ob- 
servable Universe contained within our particle horizon. 
Critically, the numbers that we consider to be "cosmic 
parameters" , for example the total matter density pa- 
rameter fini, arc subject to sample variance, since we 
estimate them by making observations within our own 
causal patch. This difference between a specific realiza- 
tion and the ensemble average is called cosmic variance. 
It is most familiar in the context of the observed values of 
the multipoles of the microwave sky. Indeed, the cosmic 
variance in those multipoles places corresponding limits 
on our ability to constrain the cosmological parameters 
using the observed multipoles. 

This decade's "era of precision cosmology" has seen 
these parameters measured to impressive precision [2]. 
Today the uncertainty in our estimates of cosmological 
parameters is dominated by the limitations of our instru- 
ments, but cosmic variance will become dominant in the 
distant future if the sensitivity of our instruments con- 
tinues to improve. In particular, experiments to map 
the distribution of matter at high redshift using 21 cm 
hydrogen line observations may be able to constrain the 
curvature parameter to within ^ 10~^ [•!]. 
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In this article, we discuss cosmic variance as it relates 
to the density and curvature parameters Om, f2i-, and 
^Ik- Ref. [4] addresses some related issues regarding the 
cosmic microwave background dipole and monopole. 



A. Overview 

The total matter content of a homogeneous and 
isotropic universe is characterized by the dimensionlcss 
time-dependent parameter 
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where H is the Hubble rate and pm is the energy den- 
sity of matter (dark matter and baryons). The radiation 
density parameter fij- is defined analogously, in terms of 
H and radiation density pr (photons and neutrinos). 
Likewise, spatial curvature is characterized by 



(2) 



where '^^^i? is the Ricci curvature scalar of the homoge- 
neous spatial slices. Thus VLk is determined by the ratio 
of the Hubble length to the spatial curvature length. The 
parameter U,k is also time-dependent: Since the spatial 
curvature decays like in an expanding universe, for 
scale factor a, we have 



nK{t)a^{t)H^{t) = const. 



(3) 



This time dependence encodes the familiar behavior that 
il.K = is unstable in an expanding matter- or radiation- 
dominated universe, leading to the flatness problem of 
hot big bang cosmologies, but stable in an inflating uni- 
verse, leading to the widely accepted resolution of that 
problem. 
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Indeed, inflation is usually thought to last much longer 
than is needed to flatten the Universe to within observed 
bounds on Uk- (Current constraints are —0.0175 < 
< 0.0085 at 95% confidence [2, Table 2].) There- 
fore it is expected that the Universe is essentially flat on 
scales much larger than the current radius of the last scat- 
tering surface, tls- Equivalcntly, the ensemble average 
(over realizations of the primordial fluctuations) of the 
spatial curvature averaged over our observable volume is 
expected to be extremely small. 

However, as it drives the spatial curvature towards 
zero, inflation also generates the near-scale-invariant 
spectrum of primordial fluctuations required to form 
structure [5] . The primordial comoving curvature pertur- 
bations carry a characteristic dimensionless amplitude of 
order 10~^, which suggests that even an ideal measure- 
ment of spatial curvature will be subject to irreducible 
random noise, or cosmic variance, of order flfc ^ 10~'^ 
(see, e.g., Knox [G]). For adiabatic fluctuations, the Ein- 
stein energy constraint equation then suggests that the 
parameters fl^ and fl^ should be subject to cosmic vari- 
ance of similar relative size. Given a primordial spec- 
trum, it is possible to use linear perturbation theory to 
quantify these expectations and calculate the cosmic vari- 
ance of, and hence the best possible constraints on, the 
density parameters flm, ^r, and fix- 

To do this, we must generalize Eqs. (1) and (2) to 
the case of universes that depart from exact homogene- 
ity and isotropy. The most natural approach is to 
take these equations to define local parameters flraix,t), 
rir(a;,t), and r2/^(a;,i) in terms of similarly local densi- 
ties pm{x,t) and pr{x,t), curvature '■■^■'i?(a;, i), and Hub- 
ble rate H{x,t),^ and to take fim, f^r and flu to be the 
spatial averages of nm{x,t), ni-{x,t), and flK{x,t), re- 
spectively, for a fixed time t. However, since we no longer 
have homogeneous hypersurfaces naturally singled out, 
we are free to vary the choice of time coordinate. For 
example, we could choose the slices of constant t to be 
spatially flat, in which case ^^^R and hence ilx trivially 
vanish everywhere! This inevitable gauge ambiguity can 
be defeated by fixing the slicing by a physical criterion. 
Since we will wish to examine the largest observable slice 
in order to obtain the best constraints, we will choose to 
work on the hypcrsurfacc of last scattering, which is de- 
fined locally by the condition that the radiation energy 
density be some constant value. 



B. Structure, assumptions, and conventions 

The structure of this paper is as follows: In section II, 
we introduce the basics of cosmological perturbation the- 



^ We define H{x, t) to be the expansion rate of the worldhnes 
comoving with matter (cold dark matter plus baryons), i.e. the 
worldlincs which observe zero total matter momentum density. 



ory, describe the problem of gauge ambiguity in defining 
cosmological parameters, and propose a solution involv- 
ing choosing new parameters defined on the hypersurface 
of last scattering. In section III, we derive the rule for 
transforming from one hypersurface to another and the 
relationship between density perturbations and density 
parameter perturbations, for matter, radiation, and cur- 
vature. In section IV, we numerically compute the cosmic 
variance of ideal measurements of the matter, radiation, 
and curvature density parameters we have defined. In 
section V, we employ an alternate method to compute 
the variance of the curvature parameter, highlighting the 
importance of the definition of curvature. In section VI, 
we adapt the results of the previous sections to give the 
variance of Qk in the present era, allowing comparison to 
observations and showing that it is impossible to measure 
^Ik to better precision than 1.5 x 10^^. 

Motivated by inflation and by current data, we assume 
vanishing background spatial curvature and strictly adia- 
batic perturbation modes. Throughout this paper we use 
the convention c = 1, and all distances are comoving. Of 
course, we also assume the Standard Model of cosmology 
to be correct. This is sensible, since we are trying to find 
minimum theoretical uncertainties, and if the real map 
of the Universe has dragons beyond our particle horizon, 
then our lower bound will certainly still hold. 



II. COSMOLOGICAL PERTURBATION 
THEORY 

The Friedmann-Lemaitre- Robertson- Walker model de- 
scribes a universe which is homogeneous and isotropic. 
Cosmological perturbation theory is the means to gener- 
alize this model to the relaxed requirement of statistical, 
rather than exact, homogeneity and isotropy. We begin 
with a brief summary of gauge choice in perturbation 
theory. See, e.g., Mukhanov [7, section 7.1] for more de- 
tails. 

We can introduce perturbations to an FLRW universe 
by writing the metric as 

Here Sg^i, is a symmetric 4-by-4 matrix, not a tensor, 
and it is convenient to decompose it as 

(4) 

breaking down the 10 degrees of freedom of into 

• four scalar degrees in the form of four scalar fields 
(f), tjj, B, and E, which give rise to structure in the 
Universe; 

• four vector degrees in the form of two divergence- 
free spatial vectors Fi and Si, which decay quickly 
as space expands; 
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• and two tensor degrees in the divergence-free, trace- 
free spatial tensor hij, which give rise to the two 
polarizations of gravitational radiation. 

The transformations of Sg^^, under coordinate changes 
are defined by the invariance of the line element ds^. In 
particular, there is a coordinate freedom, called gauge 
freedom, in mapping from the unperturbed spacetime to 
the perturbed spacetime. For example, a given perturbed 
universe can be described using coordinates in which the 
spatial hypcrsurfaccs of constant time coordinate have 
uniform total energy density, or alternatively by coordi- 
nates in which those hypersurfaces' normal curves have 
uniform local rate of expansion. The two choices result in 
different matrices dg^,y. Critical to the decomposition of 
Eq. (4) is that the three classes of perturbations (scalar, 
vector, and tensor) do not mix with each other under the 
gauge transformation from one choice of coordinates to 
another. 

Concretely, a gauge transformation in the perturbed 
spacetime can be written as a change of coordinates 

xf" ^xt" ^ x'' +^^', 

which has four degrees of freedom. So four out of the ten 
degrees of freedom of Sg^^, are unphysical gauge degrees. 
The remaining six physical degrees of freedom may be 
decomposed into 

• two scalar degrees, representing density perturba- 
tions; 

• two vector degrees, which are not of cosmological 
interest because they decay very quickly; 

• and two tensor degrees, representing gravitational 
radiation. 

Of interest to us are the scalar degrees, which may 
be expressed as linear combinations of density ("adia- 
batic" ) perturbations and entropy ( "isocurvature" ) per- 
turbations [8]. When the anisotropic stress vanishes (as 
is the case during matter domination, for example) the 
scalar sector reduces to a single physical degree of free- 
dom. 

Linear cosmological perturbation theory applies 
through most of the history of the Universe, and it ap- 
plies today on sufficiently large scales; it breaks down 
only when and where perturbations grow in magnitude 
to order unity. 



A. Power spectra and cosmic variance 

The principal connection between cosmological pertur- 
bation theory and our observations of the perturbed uni- 
verse is through the statistical power spectra of observ- 
able quantities. 

Assuming statistical homogeneity and isotropy, a 
scalar perturbation / on a given spatial hypersurface may 



be written as an isotropic Gaussian random scalar field on 
R!^ with a power spectrum P{k). That is, if the Fourier 
transform of / is given by 
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(27r)3/2 



d^k f{k)e' 



then for each fc, /(fc) is a Gaussian random variable with 
mean (over the statistical ensemble) 



/>))-0 
and variance specified by 

'f{k)r{k')) = 6^k-k')Pik), 



(5) 



depending only on k = \k\. Note that this definition 
of P{k) depends on the choice of normalization of the 
Fourier transform; the choice we have made is common 
in the literature (see Mukhanov [7, section 8.1] and Liddle 
and Lyth ["i, section 4.3.2]) and matches that used by the 
numerical tool CAMB [9]. 

P{k) has dimension k^f^, or x^f^. In particular, if / is 
dimensionless, then P{k) has dimensions of volume, and 
it is therefore convenient to define a dimensionless power 
spectrum [5, ibidem] 

V{k) ^ ^P{k). 

Critically, however, the power spectrum of / depends 
on the spatial hypersurface we have chosen to define /. If 
the hypersurface corresponds to a fixed coordinate time, 
then this choice depends on the gauge we are using. For 
example, if / is the matter density perturbation Spm/ Pm, 
then in the uniform matter density gauge, 'P(fc) = for 
ah k. 

Suppose we are interested in the variance of an es- 
timate of a dimensionless cosmic parameter F due to 
observation over a finite region V of some spacelike hy- 
persurface. To estimate the mean value F, we compute 
the mean Fy over the volume V of space in which we 
can observe F. We want the statistical variance of this 
estimate. 

Let us define a perturbation variable f = F — F . Then 
the variance of the estimate Fr (the mean over a spherical 
region of comoving radius r) is equal to the variance of 

fr- If the dimensionless power spectrum of / is 'P(fc), 
then (as we show in Appendix A) , this variance is simply 



Var 



fr 



J ^r(^k)w^^k), 



(6) 



where we have defined the window function 

ji(fcr) 



Wr{k) 



kr 
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B. Gauge and cosmic parameters 

When we describe our Universe as a perturbed FLRW 
universe, the background FLRW spacetime has weU- 
defined cosmological parameters which are spatially con- 
stant, including background values of flx{x,t) (where 
X here may stand for matter, radiation, curvature, or 
even dark energy), which we call Qx- But gauge free- 
dom tells us that there is no unique choice for the point 
identification map between the background FLRW uni- 
verse and our perturbed universe; in other words, there is 
no unique description of the real Universe as a particular 
perturbed FLRW universe. We have no preferred global 
choice of time coordinate, and since the local parameters 
fix evolve over time, wc see that perturbations from the 
background values fix, which wc will require to deter- 
mine variances, are also gauge-dependent quantities. 

It would seem as though gauge ambiguity leaves us 
with no place to plant our feet. But the real problem is 
that in a perturbed FLRW universe, these background 
parameters are simply not well-defined quantities. 

To escape gauge ambiguity, we need cosmic parame- 
ters which are independent of the time of observation 
for each FLRW universe. To retain the utility of fix, 
the simplest choice is to pick the value of Qx at some 
physically fixed time. The time we choose is that of last 
scattering, which at background level we may define as 
an instantaneous event taking place at redshift z = 1100, 
and at perturbative level as a spacelike hypersurface of 
constant radiation energy density.^ We denote the back- 
ground value of rix at last scattering by fl^- As with 
conventional parameters ^x, this definition carries over 
locally to the real perturbed Universe, although only on 
the hypersurface of last scattering itself. We may extend 
the definition in a continuous fashion to points not on the 
last scattering hypersurface, for example by defining iV^ 
at any event to be the mean of the local fix over the in- 
tersection of the hypersurface of last scattering with the 
light cone of that event, but this is not strictly necessary. 

There are other reasons that it is a good choice to 
work on the last scattering surface. At z = 1100, nonlin- 
ear structure has not yet formed, and linear cosmologi- 
cal perturbation theory is valid on virtually all distance 
scales, particularly those which dominate the observable 
perturbation power spectra. In addition, the contribu- 
tion from dark energy is negligible at this time, so we 
need consider only matter, radiation, and curvature in 
our calculations. 

fJ^if is a good choice of parameter for one more reason. 
We are interested in the theoretical limit of our ability 
to estimate cosmological parameters through electromag- 



^ Last scattering actually occurs over an interval of time as the 
baryons gradually recombine. We can define the time of last 
scattering more precisely to correspond to the maximum of the 
Thomson scattering visibility function. 



netic observations. Thus, consider an observer at event 
O who can collect complete information from their causal 
past as far back as last scattering. Let Eo be the inter- 
section of the causal past of O with the hypersurface of 
last scattering; we call this the causal last scattering hy- 
persurface of O and denote it by Eq . The region of space- 
time observable by O is the future Cauchy development 
of Eq, in which the physics is completely determined by 
information on Eq. Thus, for a given component X, the 
observer's best estimate ft^ of the background value of 
fl^ is simply the mean of fix over Eo, 

Note that, while the directly visible last scattering sur- 
face corresponds to just the boundary of Eo, wc are in- 
terested here in ultimate limits. The interior of Eo can 
be observed, for example, by mapping the distribution of 
matter with galaxy or 21 cm line observations. 

We can then compute the cosmic variance of the esti- 
mate fl^ using Eq. (6), if we know the power spectrum 
of the perturbations 5flx of fix on the hypersurface of 
last scattering. 



III. PERTURBATIONS AT LAST SCATTERING 

The evolution of a realistic set of coupled cosmologi- 
cal perturbations, the crux of cosmological perturbation 
theory, is complicated. Fortunately there exist several 
numerical tools which allow us to obtain power spec- 
tra. In this section we use CAMB [10],^ which performs 
calculations in the gauge in which the normals to the 
constant-time slices are comoving with cold dark mat- 
ter (CDM). Since CDM is assumed not to interact with 
the other matter components, this gauge is also a syn- 
chronous gauge. 

The constant-time hypersurfaces of the CDM- 
comoving gauge are not hypersurfaces of constant radi- 
ation energy density, which defines the hypersurface of 
last scattering. At linear order, the transformation of a 
four-scalar between these two slices is given by the first 
term in a Taylor expansion in proper time t. Writing X 
for scalar quantities on the CDM-comoving hypersurface 
and X for those quantities on the hypersurface of last 
scattering (note that we are now changing the meaning 
of the tilde), we have 

X{x) = X{x) + At{x)dtX{x), (8) 

where At is the proper time displacement between the 
two slices. This transformation applies to an energy den- 
sity, X = p, even though p is not a true four-scalar, since 



^ CAMB is available at http://camb.info/. 
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at linear order the density does not change under boosts. 
Similarly, it applies to the expansion rate, X = H, when 
we keep fixed the worldlines whose expansion H mea- 
sures. 

The last scattering hypersurface is a hypersurface of 
constant radiation energy density. Assuming adiabatic 
perturbations, this is also a hypersurface of constant mat- 
ter density on large (supcr-Hubblc) scales. We will see 
later that most of the contribution to the variances we 
calculate comes from scales that are super-Hubble at last 
scattering, so that it is a reasonable approximation to 
take last scattering to occur at constant matter density. 
Letting X = p^, the left hand side of Eq. (8) is constant 
in X. For the right hand side also to be constant, we 
require 



At 



dtPm{x) ■ 



where Sp^ is in CDM-comoving gauge. Taking expecta- 
tions, Eq. (8) now gives 



X -X = (AtdtX) 



(Sp..) §f ^ 0, 



where we have replaced dtX and dtPm by the mean values 
dtX and dtpy^, since we are working to first order in 
(5pm- We may therefore convert Eq. (8) into a form which 
applies to perturbations 5X = X — X: 

SX{x) = 5X{x) ^^^^dtX. 

Again replacing dtPm by dtp^, we take the Fourier trans- 
form to obtain 



6X{k) = SX{k) - ^-P01dtX. 



(9) 



In what follows, we will work exclusively in Fourier 
space. 



A. Matter density parameter perturbations SQm 

From Eq. (1), it is easy to work out a relationship 
between Sflra, Spm, and SH: 



Spn 
Pm 



SH 

) 

' H ■ 



(10) 



To calculate the perturbations to rim on the last scat- 
tering hypersurface, we first evaluate this equation on 
a CDM-comoving hypersurface and then use Eq. (9) to 
switch to a surface of constant energy density. This is 
valid since, as explained above, both pm and H can be 
treated as four-scalars. 

The first thing to do is to express SH in terms of Sp^. 
Using the continuity equation 



dtPn 



-iHp„ 



which holds to linear order when H measures the expan- 
sion rate of the comoving matter worldlines, we find 



SH = H - 
1 

^ "3 



H 

dtPn 



+ dfSpm dtp„ 



Pni 

Sp„ 

Pm 



Spn 



Pn 



dropping terms in {Spm) along the way. 

The time derivative here cannot be found analyti- 
cally in general, since the perturbation equations are 
second-order in p^. However, as we will show in Ap- 
pendix B, this information can be extracted numerically 
from CAMB. 

Finally, we employ Eq. (9) to move to the last scatter- 
ing hypersurface: 

^^m _ SQra _ S p^y dt^ra 
^in ^tPra 

Using the background relation il„i + f^r = 1 and the 
continuity equations for pm and pi-, we can rewrite the 
background component of the second term as 



1 



1 

3Pm 



Putting this all together, we have 



SVL^ 

17 rn 



1 



1- 



Spn 
Pm 



^ fS_p, 

3H V Pm 



(11) 



This equation relates the amplitude of perturbations in 
rini on the hypersurface of last scattering to the am- 
plitude of perturbations in p,,! on the nearby CDM- 
comoving hypersurface used by CAMB. 

If we make the simplifying assumption of matter dom- 
ination at last scattering, fl^ = 1, then we have the 
well-known result Spm/pm oc a (see, for example. Press 
and Vishniac [11, Eqs. 32c and 9]), so 



Pm 



= H 



Spn 



Also, the correction term for moving to the hypersurface 
of constant density vanishes. In this case, Eq. (11) is 
simply 



50n 



5 (5pn 
3 Pm 



B. Radiation density parameter perturbations 50,^ 

We can also examine the variance of fJi-, since Eq. (10) 
still holds if we substitute radiation for matter: 



sn. 



Spr 

Pr 



H 
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On the hypersurface of last scattering, adiabaticity tells 
us that 5p^ = (5pm = (on large scales), so 



The scalar index determines the departure from scale 
invariance in the primordial spectrum via 



5n, _ _ sn,. 



(12) 



(Recall that in defining both and Or, we take H to be 
the expansion rate of the same comoving matter world- 
lines.) 



C. The curvature parameter Q,k 

The spatial curvature scalar ('^^i?, for hypersurfaces or- 
thogonal to the matter comoving worldlines, enters the 
local energy constraint (Friedmann) equation as (see, 
e.g., Ref. [r,, Eqs 4.165 and 14.132]) 



(Pn 
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Thus, even with our local definitions of the parameters, 
we have the familiar relation + Oj. -I- O^ = 1. Also, 
Eq. (9) assures us that with fix = 0, O^ is independent 
of the slicing we are working on. However, it is simplest 
to continue working on the constant energy density hy- 
persurface, obtaining 



K 



1-a, 



5o„ 



(13) 



using Eq. (12). 

It is important to note that, regardless of which slic- 
ing we are using, Q.^ is defined in terms of comoving 
curvature, that is, the spatial curvature of hypersurfaces 
orthogonal to the matter-comoving worldlines. 



IV. THE VARIANCES OF , Q^^ , AND Q?^ 

In this section, we drop the use of the bar to denote 
mean values, in order to avoid cumbersome notation. 

To calculate numerical values for the variances, we 
must choose a set of background cosmological parame- 
ters. We use the following maximum-likelihood values 
from the five-year WMAP data combined with baryon 
acoustic oscillations and Type la supernovae [2, Table 1]: 

Ob,o/io = 0.02263 
^^m,o/^o = 0.1362 

/lo = 0.703 
Vn{ko) = 2.42 x 10"^ 
ris = 0.961 



Vn{k) = Vn{k^) 



k 



The corresponding value of 0^ o is 0.276. Since current 
data are consistent with (and the inflationary model pre- 
dicts) a flat 0/^ = universe, we use this value for sim- 
plicity. We also have Oj-^o ~ 4.17 x lQ~^h^^ [■■<, section 
2.2.2]. Together, these give Oa,o = 0.724. 

With these parameter values, the comoving distance to 
last scattering is tls ~ 9.86/i~"'^ Gpc ~ 14.0 Gpc, and at 
last scattering we have « 0.252 and Oj^f sa 0.748, so 
it is not adequate to assume matter domination at last 
scattering (see Appendix C for the details). 

Using Eqs. (6) and (11) together with data produced 
using CAMB with the above parameter values as input 
(see Appendix B for the details), we have computed the 
fractional variance of Vlltf to be 



Var ( 1 = 1.51 X 10-^^ 



The absolute variance of Ojl^ is therefore 



Var 



= in]'^) Var 



= 8.44 X 10' 



Its standard deviation, which we shall call , is 



LS 



Var OLS) =2.90x 10" 



By Eq. (12), the fractional variance of Q}^^ is equal to 
that of ^llif : its standard deviation is thus 



Var {0\ = 9.79 x 10^. 



Finally, using Eq. (13), we see that the absolute vari 
ance of is equal to the fractional variance of Oj^f 
its standard deviation is 



so 



LS 



X 10" 



(14) 



Note that a^? + a. 



LS 



(7^ , since perturbations in 



matter and radiation are perfectly correlated by the as- 
sumption of adiabaticity. 



DIRECT CALCULATION FOR ^a- 



Here subscript refers to values today, h is the Hubble 
rate in units of lOOkms"^ Mpc^^, VvSJ^) is the power 
spectrum of the primordial comoving curvature perturba- 
tion Tl, and WMAP uses a pivot scale fco = 0.002 Mpc"^ 



In the preceding sections we arrived at a result for the 
variance of VlS^ via a calculation for Om and the Einstein 
energy constraint equation. It is possible to take an al- 
ternative and more direct route to , working straight 
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from the definition Eq. (2). The idea is simply to calcu- 
late the spatial curvature of the last scattering hypcrsur- 
face itself, i.e. the curvature of a hypersurface of uniform 
radiation density, then average over the (in principle) ob- 
servable volume So using Eq. (7), and finally take the 
variance. 

The spatial Ricci curvature ^'^^R of the last scattering 
hypersurface is related to the metric curvature perturba- 
tion tjjry (the isotropic scalar part of the perturbed spatial 
metric) by [5, Eq. 14.129] 



(15) 



significant contribution from Hubble-scalc modes, which 
largely accounts for the discrepancy. Using COSMICS, 
we can repeat the above calculation, but computing the 
curvature of the comoving slice at last scattering instead 
of the uniform radiation density slice. The result is 



3.74 X 10" 



(20) 



which is now very close to Eq. (14). We believe the 
remaining disparity is due to the fact that the CDM- 
comoving gauge used by CAMB does not coincide pre- 
cisely with the total matter comoving gauge. 



where ip^ is the scalar -0 of Eq. (4), evaluated in the uni- 
form radiation density gauge (here radiation refers only 
to photons, not neutrinos, hence the subscript 7). Since 
we wish to write the variance in terms of the power spec- 
trum of the primordial comoving curvature perturbation 
TZ, we define a transfer function T^lk) through 



V'^(fc,tLs) = ^x(fc)7^(fc). 

Combining these expressions with Eq. (2) gives 

2 



(16) 



3 V^LS-ffLS 



TK{k)n{k), (17) 



for the curvature parameter at the time of last scatter- 
ing. Finally, averaging over the sphere So said taking 
the variance using Eq. (6), we find 



(18) 

Now all that is required is to determine the transfer 
function Tft-(fc) for the uniform radiation gauge curvature 
perturbation. Such detailed information can be readily 
extracted from the software package COSMICS [12].'^ 
Extracting the transfer function and performing the in- 
tegral in Eq. (18), we find standard deviation 



LS 



4.35 X 10" 



(19) 



To understand the apparent disparity between this re- 
sult and that of Eq. (14) above, recall from section IIIC 
that the curvature parameter calculated there describes 
the spatial curvature of hypersurfaces orthogonal to 
the comoving worldlines. These comoving slices coin- 
cide with uniform total energy slices on super-Hubble 
scales [13], and hence with uniform radiation slices for 
adiabatic modes. The integrand in Eq. (18) receives most 
of its support from scales between the Hubble scale at last 
scattering and scales about 100 times larger, so we ex- 
pect the two estimates to be close. However, there is a 



VI. INTERPRETATION AND DISCUSSION 

Our calculations in the preceding sections are all eval- 
uated at the time of last scattering. Thus they are not 
directly comparable with standard determinations of the 
density parameters r^m^Oi ^r,Oi and ^Kfli which are eval- 
uated today. For a homogeneous background spatial cur- 
vature, it is simple to relate the curvature parameter fix 
at different times using Eq. (3). Namely, Qk is propor- 
tional to the square of the comoving Hubble length, since 
the comoving spatial curvature length remains constant. 
For the inhomogeneous case, although variances are per- 
fectly well-defined on the hypersurface of last scattering, 
we are again confronted with ambiguity in how to trans- 
late them to values today. We choose to translate the 
standard deviation of fix, Eq. (14), in the same way as 
a background curvature, using Eq. (3) to obtain 



''^ \ aoHo J 



1.58 X 10 



-5 



(21) 



where we have used the Friedmann equation to determine 

We can similarly translate the results from section V 
to today. For the spatial curvature of the last scatter- 
ing hypersurface itself, i.e. the uniform radiation density 
shce, Eq. (19), we find 



= 1.77 X 10-^ 



(22) 



Translating the curvature of comoving slices, Eq. (20), to 
today gives 



'^K.n = 1-52 X 10 



-5 



(23) 



which again is very close to the estimate Eq. (21). 

To help understand this evolution, it may be helpful 
to consider what happens to the variance of fix at fu- 
ture times, again using the prescription of Eq. (3). The 
variance will approach zero as we approach the late time 
de Sitter stage, ^ when the comoving Hubble length de- 



COSMICS is available at littp://web. mit.edu/edbert/. 



^ Here we are assuming for simplicity that the dark energy is a 
pure cosmological constant. 
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FIG. 1: The time dependence of ctk. The change of slope 
is due to the transitions from radiation domination to matter 
domination and from matter domination to dark energy dom- 
ination. The standard deviation today has already begun to 
decrease as we approach de Sitter. 



creases indefinitely; see Fig. 1. This is a refiection of pre- 
cisely the same mechanism by which inflation drives spa- 
tial curvature towards zero in the early universe. Indeed, 
the variance calculated for today in Eqs. (21) to (23) 
has already begun to decrease significantly, as we are al- 
ready approaching A-domination. However, it could be 
argued that this decrease in the variance of flu at late 
times is deceptive, in that we are simply scaling essen- 
tially the same last-scattering variance to a smaller and 
smaller comoving Hubble length, (r^s will increase some- 
what in the future, but it will not grow without bound.) 
Recall that the large-scale comoving curvature pertur- 
bations are described by a constant amplitude of order 
10~^. An alternative prescription to Eq. (3) for evolv- 
ing the curvature variance would be to scale it with tls 
rather than the Hubble length. Then the variance in il.K 
would approach a constant at late times, of order 10~^. 

In summary, we have quantified the common lore that, 
within the standard inflationary paradigm, we can never 
measure fix to better precision than of order 10"^. Our 
calculations are completely independent of the details of 
any particular observational technique and hence provide 
a fundamental limit for any measurement of curvature. 
In particular, it will be impossible to determine the sign 
of , and thus whether the Universe appears to be open 
or closed, if |Sl/i-,o| ^ cr^ ^ 10~^. On the other hand, if 
we some day determine |riA',o| 3* <^k- then we will know 
that (within the framework of the Standard Model) we 
have measured the very large scale geometry of our Uni- 
verse and not merely a local fluctuation. While current 
observational constraints on curvature are well above this 
fundamental limit, the mapping of a significant portion 
of our past light cone, as would take place with planned 
21 cm hydrogen line experiments [:>] , will bring us close to 
our ultimate ability to constrain curvature, matter con- 



tent, and other cosmological parameters. 
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APPENDIX A: SCALAR FIELD VARIANCE 

Consider the mean value of a function / : i?"^ R over 
some spherical volume ^T^r^ ( "top-hat" ) of flat space. We 
will denote this quantity by /,.: 

fr^T^ f f{x) 

|7rr3 7|^|<^ (27r)J/^ J 

- (2^)3/2 J ^ ^ nfe) ■ 

Here j\ is the spherical Bessel function of first order, 
, , sin X — X cos X 



The variance of jr is 




If / is a perturbation, then \Jrj = {^fj = 0, and thus 
the variance of fr is simply 




where we have used Eq. (5). 
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APPENDIX B: CAMB DETAILS 

Here we provide some details of our numerical imple- 
mentation for the benefit of readers who wish to repro- 
duce the results for different models. 

CAMB produces two types of data files. 

The matterpower files are two-column files in which 
the first column gives values of k/h^ and the second col- 
umn gives corresponding values of h^P{k)^ where is 
the dimcnsionlcss Hubble parameter evaluated today and 
P{k) is the dimensional power spectrum of the quantity 
Spm/Pm in the CDM-comoving gauge. 

The transfer files have seven columns, in which the 
first column is again k/ho and the remaining columns 
contain internally-defined dimensional transfer functions. 
We are interested in the total matter transfer function, 
which is found in the seventh column. We will call it 
Tc(k) to distinguish it from the dimcnsionlcss transfer 
function Tik) commonly used in the literature. 

Tc is related to the matter power spectrum P{k) by 

P{k)^T^{k)-2^-'kVn{k), 

where Vnik) is the dimcnsionlcss power spectrum of the 
primordial comoving curvature perturbation, which can 
be constrained observationally. (For a scale-invariant pri- 
mordial spectrum, V-jiik) = const.) Thus wc have the 
convenient relation 



Pik) ^ k^T^{k)Vn{k). 
From this, we see that 

^ oc fc2Tc(fc), 

Pm 



(Bl) 



SO 



Q |'^Pm^ _ dtTcik) 6p„ 




Pm ) Tcik) 
We may therefore rewrite Eq. (11) as 



FIG. 2: The integrand components VF^^^s ^^'^ T^k^TcVn- 
Note the rapid oscillations of W relative to the spacing of 
the data points. 



To evaluate this expression numerically, it is not suffi- 
cient simply to perform a Riemann sum, since the win- 
dow function IFrLs(^) oscillates on a scale of r^g ~ 
10-''Mpc"\ but CAMB data points can be separated 
by up to ~ 10-2 Mpc"^ (The function Tc{k) varies 
much more slowly with k than does Wrj^g{k).) The so- 
lution is to construct linear interpolations of F(fc)Tc(fc) 
and then to perform a piecewise integral, including tail 
terms. We have done this on a set of data produced by 
CAMB using the parameter values given in section IV 
and with the CAMB parameters transfer Jonax set to 
1 and accuracy_boost set to 3. 

We observe that typical values of F(fc) are close to 
2, representing a correction of 20% to the matter- 
domination value of 5/3. The integrand components 
and T^k'^T^Vn arc plotted in Fig. 2. 



APPENDIX C: PARAMETERS FOR LAST 
SCATTERING 



Sii=F(fc)^ 



(B2) 



To estimate fir at last scattering, we simply work back- 
ward from today: 



where we have defined 

r(fc).i + la. + ^^. 

^ ^ 3 ' 3H Tc{k) 

Hence, we have modified the CAMB source code to out- 
put drtt and drTcik) as well (where t is conformal time), 
allowing us to compute T{k) numerically. 

Inserting Eqs. (Bl) and (B2) into Eq. (6), we see that 
the fractional variance of is 



dk 



T\k)k^T^{k)Pnik)W^,Ak). 



Pr,oaLS 



I 1 , ^rn,0 
= I 1 + — «LS 
"r,0 



^A.O 4 
"r,0 
W 0.252, 

with ao = 1 today. Similarly, we find 
n}:^ « 0.748, 



and f^A at last scattering is negligible. 

To calculate the comoving distance to last scattering. 
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wG start with the form of the Friedmann equation 

(a) " (^r.oa^"^ + ^m,oa-'^ + ^a,q) 

in terms of density parameters evaluated today, with the 
dot denoting differentiation with respect to proper time. 
Therefore the distance to last scattering is 

/■*" dt _ n' da 
_ 1 /"^ da 

Hq + Om,oa + (1 - f^r,0 - ^m,o) 

« 9Mh-^ Gpc. 
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